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Abstract 

We present the next-to-leading order partonic cross sections involving an elec- 
troweak vector boson and four massless partons (quarks, gluons or long living 
gluinos) in a general gauge theory with a simple Lie Group. The vector boson couples 
to a massless lepton pair and a quark-antiquark pair. The cross sections are given 
in terms of group independent kinematical functions multiplying the eigenvalues of 
the Casimir operators of the Lie group. This kind of color decomposition is required 
for the calculation of 0(aj?) corrections to the group independent kinematical func- 
tions in the four-jet production cross sections in electron-positron annihilation. The 
knowledge of these corrections facilitates the simultaneous precision meaurement 
of the strong coupling and the color charge factors using the four-jet LEP or SLC 
data as well as the test whether these data favour or exclude the existence of a light 
gluino. 

PACS: 12.38Bx, 13.38.Dg, 11.30.Pb 



Electron-positron annihilation into four jets is a very clean way of measur- 
ing the basic properties of Quantum Chromodynamics (QCD), the theory of 
strong interactions. This process allows the simultaneous measurement of the 
strong coupling a s and the eigenvalues of the Casimir operators of the under- 
lying symmetry group, the color factors. Such a measurement provides a very 
stringent test of QCD because a s is the only free parameter in the theory, 
while the color factors show whether the dynamics is indeed described by an 
SU(3) symmetry. However, in order to carry out such a full measurement the 
theoretical prediction of perturbative QCD has to be known at next-to-leading 
order accuracy. Lacking this sort of precision of the theoretical description, so 
far the experiments used four-jet data for measuring the ratios of the color 
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factors Ca/Cf and Tr/Cf that are relatively insensitive to higher order correc- 
tions [1]. The a s measurements require the fixing the normalization of four-jet 
observables, which however, suffers large radiative corrections. For example, 
the experimental four-jet rate measured as a function of the jet resolution 
parameter was about a factor of two higher than the leading order prediction. 

Recent theoretical developments make possible the next-to-leading order cal- 
culation of four-jet quantities. There are now several general methods available 
for the cancellation of infrared divergences that can be used for setting up a 
Monte Carlo evaluation of next-to- leading order partonic cross sections [2]. 
The other vital piece of new information is the one-loop matrix elements for 
the relevant QCD subprocesses — i.e. for the production of four-quarks and for 
the production of two quarks and two gluons in electron-positron annihilation 
- that are now available [3-6]. In refs. [3,5] Campbell, Glover and Miller make 
FORTRAN programs for the next-to-leading order squared matrix elements 
of the e + e~ — > 7* — > qqQQ and qqgg processes publicly available. In refs. [4,6] 
Bern, Dixon, Kosower and Wienzierl give analytic formula? for the helicity 
amplitudes of the same processes with the e + e~ — > Z° — > four partons chan- 
nel included as well. The helicity amplitudes for the processes e + e~ — > qqggg 
and e + e~ — > qqQQg have been known for a long time [7]. Using the helicity 
amplitudes in refs. [4,6,7], Dixon and Signer calculated the next-to-leading 
order corrections to four-jet fractions for various clustering algorithms [8], as 
well as to angular distributions [9] using a general purpose partonic Monte 
Carlo program [10]. As a result, the agreement between the data and the 
theoretical prediction for the four-jet rates improved significantly. In the cal- 
culation of these four-jet fractions the underlying gauge group is assumed to 
be SXJ(Nc) and the results are presented in the form of an expansion in the 
number of colors. This way of presenting the results is sufficient if one is inter- 
ested in the four-jet rate and will certainly find important application in LEP2 
physics, where the e + e~ — > (7, Z) — > 4 jets events are the main background 
to e + e~ — > W + W~ — > 4 jets events, but not for increasing the precision in 
measuring the color charge factors. 

One can measure the color charges directly using any kind of four-jet ob- 
servable, provided the partonic cross section is given in a factorized form of 
eigenvalues of the Casimir operators and group indpendent kinematical fac- 
tors. Such a decomposition has been known for a long time at tree level [11]. 
In this letter we give the necessary decomposition of the virtual corrections. 
Combining these with similarly decomposed real corrections, we can obtain 
the full O(a^) corrections to the group independent kinematical functions in 
the four-jet production cross sections in electron-positron annihilation [12]. 

An additional piece of motivation for using a group independent color decom- 
position comes from the suggestion that the measurement of the color charges 
may be influenced by the production of light, colored, but electrically neutral 
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fermionic particles, such as gluinos [13]. The effect of the presence of light 
gluinos on the strong coupling is to increase a s (Mz) by 2% [14]. It would 
be interesting to see if the presence of light gluinos influence the value of the 
color charges more noticably. In order to make possible a next-to-leading or- 
der analysis of this effect, we also present the one- loop helicity amplitudes for 
the e + e~ — > (7, Z) — > qqgg process. We have used these amplitudes to show 
that in order to obtain the most significant light gluino exclusion limit from 
four-jet data, it is advantageous to define jets in the range of small resolution 
parameter [15]. For example, in an analysis based upon ALEPH data for four 
jet rates [16] using Durham clustering algorithm [17] at y cut = 0.002 we were 
able to exclude the existence of light gluinos at the 95% confidence level [15]. 

In presenting the new color decomposition, we rely on the work of Bern, Dixon, 
Kosower and Wienzierl who published the explicit form of the one-loop helicity 
amplitudes for the QCD subprocesses [4,6]. We use the same notation as these 
articles and introduce new ones but to the extent that is necessary. 

Our aim is to give the next-to-leading order squared matrix elements for the 
e + e~ — >■ qqpp processes (p = Q, g, or g) in terms of color factors multiplied by 
group independent kinematic functions. In order to find a group independent 
decomposition of the squared matrix element, we have to give different color 
decompositions of the one-loop helicity amplitudes for the various processes 
than those presented in refs. [4,6], where the color charge information has 
been lost by the use of the SU(iVc) Fierz identity and SV(Nc) relations Cf = 
(Nc 2 — 1)/N C , Ca = 2N C (for T R = 1). In the new color decomposition we 
can only use the defining relation of the Lie algebra 

N A 

[T a ,T% J = iY / f abc T^, (1) 

c=l 

and the definition of the quadratic Casimir invariants Cf, Ca and T R : 

Na N a 

J2(T a T^% = C F 5 tJ , J2( FaFla U = C A S cd , Tr(T«Tt b ) = T R 5 ab ,(2) 

a=l a=l 

where Na is the dimension of the adjoint representation of the gauge group, 
T a are the generators in the fundamental and F a are those in the adjoint rep- 
resentation. The latter are related to the structure constans via F£ c = — if a b c - 
Careful analysis of the color structure of the individual Feynman diagrams 
contributing to a given process shows that the color charge information can 
completely be recovered from the primitive amplitudes of refs. [4,6] for the 
QCD subprocesses, while in the case of two-quark two-gluino production mi- 
nor modification of the partial amplitudes is necessary. In order we can use 
those primitive amplitudes we use T R = 1 normalization. 
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Let us first consider the e+e — >■ qqQQ process. At tree level the new color 
decomposition of the *4.g ree (l g , 2q, 3q, 4 5 ; 5 g , 6 e ) helicity amplitude is 



N A 



AT e (l q , 2q, 3 q , 4,; 5 g , 6.) = 2e 2 <? 2 £ T^Tqq 

c=l 

"cJ* he i6 ; o(l h *, 2 h «, 3 h «, 4 h «; 5 fe , 6^) 



(3) 



x 



+C3 Q/le i 6; o(3 ft£ 3, 4 h «, 2 A 0; 5 he , 6 



ho /I ho 1 ho ohn.Khe eh. 



where e is the QED coupling, g the QCD coupling h x is the helicity of particle 
x and the coupling factors C^ he depend on the charge Q Q of quark Q in units 
of e and on the left- and right-handed couplings of the electron {y e L R ) and of 
the quark (v2 R ) as follows: 



C++ = -QC + v%v R V z (s 56 ) , C+- 

c Q + = -g s + ^Pz(s 5 6), c Q - 



-g s + ^2P z ( S56 ), (4) 
-g s + ^£p z ( S56 ). (5) 



The left- and right-handed couplings of the electron and quark are the stan- 
dard ones as given in ref. [4]. Aq-q denotes the tree-level partial amplitude 
Ag ree of ref. [4], where it was defined to include a photon propagator. The 
ratio Vz(s$q) = 556/(556 — M% + iT z M z ) appearing in the coupling factors 
replaces this photon propagator with a Z propagator [Mz and Y z are the 
mass and width of the Z° boson). 

At one loop the new color decomposition of the helicity amplitudes is 



^ lo ° P (l g ,2Q,3 Q ,4 9 -;5 e -,6 e ) = 2e^ 



2„4 



x IC. 



h„h e 



N A 



i 6;1 (l\ 2*«, 3*« 4**; 5" g , 6 he ) £ T ^ T QQ 

c=l 
N A 

+A 6 , 2 (l h \2 h Q, 3 h Q,A hq ; 5 he ,6 he ) £ (T c T d ) qq (T d T c ) ( 



+C, 



h Q h e 



c,d=l 
N A 



i 6;1 (3^, 4 h % l h ",2 h ; 5 h % 6^) £ T^Tqq 



c=l 
Na 



+A 6 ; 2 (3 hQ A hq A hq ,2 hQ ;5 hg ,6 he ) E {T c T d ) q q(T d T c ) l 



c,d=l 



N A 



(6) 



+C£i 6;3 (l'\ 2 h Q, 3 h «, 4"«; 5 he , 6 Ae ) ]T , 
where vanishes for the photon and W ± boson, while for the e+e~ — > Z° 
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qqQQ process it is 

with 9 W being the Weinberg angle. In eq. (6) we used the notation A 6;i for 
the partial amplitudes in order to distinguish them from the functions 
introduced in ref. [4], where the basic gauge invariant classes of colorless am- 
plitudes, the 'primitive amplitudes' are also given explicitly. Our new partial 
amplitudes can be expressed in terms of the same primitive amplitudes multi- 
plied by color factors. The expressions depend on the helicities of the partons. 
There are only two independent helicity configurations, which can be taken to 
be (1+ , 2q, 3q, 4^; 5g , 6+), and from which the amplitudes of the other helicity 
configurations can be obtained [4]. 

Formulas (3) and (6) apply to the case of unequal quark flavors, q ^ Q. The 
equal flavor amplitude may be obtained from the unequal-flavor formula by 
subtracting the same formula with q and Q exchanged and then setting Q = q 
in all the coupling constant prefactors [4]. 

The explicit expressions for the A 6;i partial amplitudes in terms of primitive 
amplitudes are (we suppress the 5, 6 labels of the lepton pair) 

i 6;1 (l+,2j,3Q,4r)= (8) 

-C F A*(2, 3, 1,4) + ^ (4(2, 3, 1, 4) - Af (1, 3, 2, 4)) 

+4< ++ (l, 2, 3, 4) - N f 2, 3, 4) + 2, 3, 4)) , 

i 6;2 (l+ 2J, 3q, 47) = A+ + (l, 2, 3, 4) + A+~(l, 3, 2, 4) , (9) 
i 6;3 (l+2j,3 Q ,4r)=^ x (l,4,2,3), (10) 

while for the other helicity configuration 

i 6;1 (l+,2^,3+,4r) = (11) 

C F Af(3, 2, 1,4)-^ (4(3, 2, 1, 4) + A+ + (l, 3, 2, 4)) 
+4 + ~(l, 2, 3, 4) - N f (4 + "(l, 2, 3, 4) + ^+"(1, 2, 3, 4)) , 
i 6;2 (l+ 23, 3+ 4r) = 3, 2, 4) + 2, 3, 4) , (12) 

i 6;3 (l+ 2 , 3J, 4r) = -^-(1, 4, 3, 2) . (13) 

The gluino is a majorana fermion in the adjoint representation of the gauge 
group and does not couple to the vector bosons directly. Therefore, the e + e~ — > 
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qqgg subprocesses has similar color decomposition to that of the e + e — >• qqQQ 
subprocess. At tree level this decomposition reads 



^6 3^, 4 g ; 5 e , 6 e ) — 



c=l 



2eV C^i 6;0 (l\ 2 ft2 , 3 ft3 , 4 fc «; 5 h % 6"°) ]T T^F| 2 - 3 , 
while at one-loop it is 



(14) 



Al- loop (l q ,2- g ,3 § ,4g;5 s ,6 e ) = 2e 2 g 



2„4 



x ^C. 



ho he 



N A 



i 6;1 (l\ 2» 2 , 3 h3 , 4*; 5* 6 h <) ]T 



(15) 



+ Ch 



c=l 

+i 6;2 (l^,2 h2 ,3 h3 ,4^;5 h %6^) £ (T c T d ) qq -(F d F% 3h 

c,d=l 
N A 

i 6;3 (l\ 2 h \ 3 h:i ,4 h *; 5 h % 6 fte ) ]T T c qq F~ 93 ~ g 



c=l 



The v4 6; j partial amplitudes for the e + e~ — >• qqgg process are closely related 
to those of the e + e~ — > qqQQ process. In fact, the A 6;2 and A 6;3 amplitudes 
are exactly the same, 



i 6;2 (l+, 2± 3? , 47) = i 6;2 (l+ 2g, 3S, 47) , 



9 19191 q 

i 6;3 (l+ 2± 3J, 47) = A 6;3 (l+ 2±, 35, 4") , 



> Z V"<? ' Qi "<? 
5 ' 2 0' 3 <3' 4 <? 



(16) 
(17) 



while the A 6;1 amplitudes differ in terms arising from fermion bubble and 
parent triangle diagrams [4], 



i 6;1 (l+, 2± 3J, 47) = i 6;1 (l+ 2j, 3^, 4, ) - (C F - CU)^' +± (1, 2, 3, 4) 

-N~ g ^ (4* +± (l, 2, 3, 4) + ^+±(1, 2, 3, 4)) , (18 

where N g is the number of light gluino flavors and 



A 6 A ' +± (1,2,3,4) 



c r 4 rcc ' +± (l,2,3,4) 



£ 2 \S23J 2e \-s23 



2 \ £ 



(19) 
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with cr and Aq CC ' ++ given in ref. [4]. We remark here that gluinos are Majorana 
fermions therefore, the cross section requires an identical-particle factor of | 
in the phase-space measure. 

Finally, the new color decomposition of the helicity amplitudes for the e + e~ — > 
qqgg process at tree level is given by 



2„2 



^ rcc (l 9) 2„3 9 ,4 g -;5 e -,6 e ) = 2e^ 



h„h e 



A 6 . (l hq , 2 h \ 3 hs , 4 h «; 5 h % Q he )(T 92 T 93 ) qq 
+A 6 . (l h \ 3 h:i , 2 h2 ,4 h «; 5 h % 6 he )(T 93 T 92 ) qq 



and at one-loop it is 



(20) 



^ lo ° P (l g ,2„3 s ,4 (7 ;5 e -,6 e ) = 2eV 

x [C^ he [A 6 . 1 (l h \2 h2 ,3 h3 ,4 h «-,5 h %6 he )(T 92 T 93 ) qq 
+A 6 . 2 (l h ", 2 h2 , 3 h3 ,4 h «; $ he , 6 he ){T 93 T 92 ) qq - 

N A 

+^ ; 3(l\2^3^4^5^6^) E (T c T d ) qq (F d F c ) gsg2 

c,d=l 



(21) 



+ E \ ( c l he + c v e ) ^- A {i h \ 2 h2 ,3 h3 ,A h '; $ h % Q he ) 

N A 

x E |Tr (T 92 T 93 T c )T qq + Ti(T 93 T 92 T°)T qt 



c=l 



i£ 5 (l\ 2 h2 , 3 h3 ,4 h <; 5 h % Q K ){T 92 T 93 ) qq 

+A%. 5 (l h \ 2 h2 , 3 h3 ,4 h '; 5 he , 6 he )(T 93 T 92 ) qq 
+A%. 5 (l h «, 2 h2 , 3 h3 ,4 h «; 5 h % 6 he ) 

N A 

E [Tt(T 92 T 93 T c )T^ + Tt{T 93 T 92 T c )T c 



c=l 



The partial amplitudes A Q]i can easily be constructed from the primitive am- 
plitudes of ref. [6] : 



i 6;1 (l+, 2 h2 3 h3 4~ q ) = C F A 6 (l q , 4 q ; 3, 2) 



(22) 



+^ [A 6 (l q , 2, 3, 4 q ) - A 6 (l q , 4 q ; 3, 2) - A 6 . 3 (l q , 4- q] 3, 2)] 

+ K' ++ (1„ 2, 3, 4 q ) - N f (4' ++ (l q , 2, 3, 4 q ) + A s 6 ' ++ (l q , 2, 3, 4 q )) , 
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i 6;2 (l+ 2 h g \ 3j 3 , 4r) = C^(l„ 4 ? ; 2, 3) 
+^[A 6 (l„3,2,4,-)-^(l„4,-;2,3)] 



(23) 



+A^ + (l q , 3, 2, 4 5 ) - N f (4> ++ (l q , 3, 2, 4 ? ) + A s 6 ' ++ (l q , 3, 2, 4 ? -, 
i 6;3 (l+ 2j",3^,4-) = A 6;3 (l+ 4r;2j 2 ,3j 3 ) , (24) 
A;4(l g + ,2^,3j3,4r) = ^(1+ 4-;2j",3j») , (25) 



aA (-1+ <yh 2 oh 3 A- 



a ax /i+ /i-. 9/12 q/i3\ _ 4ax/i+ A-.ohs nh 2 \ 
/± &;A\ L q i^q 1 1 °g J I± Q;A\ L q i ^q 1 ° g 1 ^ g ) 

, a ax /-i + ,(-.9/12 o/»3 A 



(26) 



IS /i+ q/i 2 a—\ — \ A ax {~\+ A-.ohs nh 2 \ /jax/i+ 4-. 9/12 0^3 ^ 

- A 6;5V ± g > Z <? > °s > % J — 9 /1 6;4V 1 g i^q^g i Z g I /1 6;4V 1 g ' ' Z ff ' °s J 



(27) 



1 r 



^(1+2^,3^4,-) = - ^ 4 (1+ 4-;2^,3^) +^ 4 (1+ 4-;3j«,2^) 



/i ax /1 + /i-. 9/12 q/i3 N 



(28) 



Before closing the discussion of the helicity amplitudes, we should remark that 
in order to take into account the effect of light gluinos in a next-to-leading 
order calculation fully, one has to make the change 



N f 



N f + 



C A 



N~ n 



(29) 



in eqs. (8), (11) and (22), (23), where in SU(iV c ) theory C A = 2N C for T R = 1. 

As a next step we use eqs. (3), (6), (14), (15), and (20), (21) to give a gauge 
independent decomposition of the next-to-leading order squared matrix ele- 
ments. The general expression for these matrix elements is 



d(T NLO, virtual = 2 J- £ Re[X ree *^6 _I ° OP ] 
helicities colors 



(30) 



We want to evaluate the color sum in such a way that the group independent 
information is maintained. For this purpose, we use the commutation relation 
(1) together with the definition of the quadratic Casimirs, eq. (2) to derive 
the necessary Lie-algebra relations. In the case of production of two unequal 
flavor quark pairs these are: 



N A 



Tr(T a T t6 )Tr(T a T tf> ) = N C C F , 



(31) 



a,6=l 
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N A 

Tr(T a T 6 T tc )Tr(T tc T 6 T a ) = C 3 , (32) 

a,6,c=l 

where in eq. (32) we introduced C3, the square of a cubic Casimir, that cannot 
be reduced to an expression of the quadratic Casimirs C F and Ca- In SU(Nc) 
C 3 takes the value (N c 2 - \){N C 2 - 2)/N c . In the case of production of two 
equal flavor quarks, in addition to relations (31) and (32), we need two more 
color sums: 



]T Tr (T a T^ b T a T^ b ) = N C C F (C F -—) , (33) 



a, 6=1 

N A / n \ 2 

(34) 



2 

Tr(T a T b T tc T b T a T tc ) = N C C F ( C F - — ) 

a,b,c=l ^ 2 ' 



When a quark and a gluino pair appears in the final state, we use relations 



N A 

J2 Tr(T a T tfe )Tr(F a F tb ) = N C C F C A , (35) 

a, 6=1 

N A /C \ 2 

Tr(T a T b T tc )Tr(F tc F fe F a ) = N C C F ( -y ) , (36) 



a,b,c=l 



Finally, for the production of a quark pair and two gluons the additional 
necessary color sums are 



N A 

J2 Tr(T a T 6 T t6 T ta ) = N C C 2 F , (37) 

a, 6=1 

N A 

(38) 



£ Tr(T a T fe T^T^)(F b F a ) S3S2 = iV c C F (^Y , 

a, 6,92,93=1 



and 



^ Tr(T a T fe T t92 T t93 )(F 6 F a ) 9392 = , (39) 

a,&,92,93=l 



2 Tr(T S2 T 93 T c )Tr(T t92 T t93 T c ) = C 3 - C F —N C . (40) 

92,93, c=l 



Using eqs. (31-40), one can evaluate the color sum in eq. (30) easily. The 
resulting differential cross sections can be written as a quadratic form: 
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— da, 
0o 



NLO, virtual 



(?) 



2tt 



(41) 



x [C (p) + C x (p) x + C„(p) y + C z (p) z 

+ C XX (P) X 2 +C X y(p)Xy + Cyy(p) ?/ 2 



In this equation o"o denotes the Born cross section for the process e + e~ — ► qq, 
p is the collection of the final state momenta, and x, y and z are ratios of 
eigenvalues of the Casimir operators defined as 



x 



Ca 
C F 



Tr 

c F 



Nc 
N A 



c 3 

N c Cl 



(42) 



These ratios are the sole quantities together with the overall normalization 
that carry group information. The coefficients Cj(p) are the group indepen- 
dent kinematical functions that depend on the subprocess. Their explicit ex- 
pressions are quite complicated, but can be obtained straightforwardly using 
the formulas of this letter. These expressions will be published in the form of 
C++ code in the partonic Monte Carlo program DEBRECEN. 

In this letter we have shown explicitly how to rewrite the one-loop helicity am- 
plitudes of the e + e~ — > (7, Z) — > 4 partons processes in a form from which the 
group independent color decomposition of the next-to-leading order squared 
matrix elements for these processes can be obtained. These kinematical func- 
tions for the QCD subprocesses can completely be written in terms of the 
primitive amplitudes given in refs. [4] and [6]. In the case of the e + e~ — > qqgg 
subprocess one has to modify the primitive amplitudes of the four-quark sub- 
process slightly. We also presented the general structure of the next-to-leading 
order partonic cross sections in terms of group independent kinematical func- 
tions multiplying ratios of eigenvalues of the Casimir operators of the Lie 
group. This kind of color decomposition together with a similar decomposi- 
tion of the e + e~ — > (7, Z) — > 5 partons squared matrix elements is required 
for the calculation of 0{a%) corrections to the group independent kinematical 
functions in the four-jet production cross sections in electron-positron annihi- 
lation. We anticipate that these corrections will improve our knowledge of the 
basic parameters of QCD. 
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